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We study antiferromagnetic spin-1/2 Heisenberg ladders, comprised of n c chains (2 < n c < 6) with 
ratio J±/J\\ of inter- to intra-chain couplings. From measurements of the correlation function we 
deduce the correlation length f(T). For even n c , the static structure factor exhibits a peak at a 
temperature below the corresponding spin gap. Results for isotropically coupled ladders (Jj_/ Jn = 1) 
are compared to those for the single chain and the square lattice. For J±/J\\ < 0.5, the correlation 
function of the two-chain ladder is in excellent agreement with analytic results from conformal field 
theory, and £(T) exhibits simple scaling behavior. 



PACS numbers: 75.10.Jm, 75.40. Mg 

Low-dimensional quantum Heisenberg antiferromag- 
nets (QHA) exhibit many fascinating properties. In 1931, 
Bethe demonstrated for the one-dimensional (ID) spin 
5 = 1/2 nearest-neighbor Heisenberg chain that quan- 
tum fluctuations prevent the existence of an ordered 
ground state [EJ. Instead, this system exhibits power-law 
correlations with gapless excitations. Haldane suggested 
in 1983 that all non-integer-spin chains are gapless, but 
that integer-spin chains should have a spin gap pi. By 
now, there is much evidence for the correctness of this 
famous conjecture |J. For the two-dimensional (2D) 
analog of the spin chains, the square-lattice nearest- 
neighbor QHA, it has been established over the past 
decade that an ordered ground state exists even in the 
extreme quantum limit of S = 1/2 . 

Spin ladders are arrays of coupled chains, and thus 
present interpolating structures between ID and 2D ||. 
These systems are thought to be realized, with 5 = 1/2, 
in the materials (VO)2P20j || and Sr n -\Cu n+ i02n @- 
They allow the study of the dimensional crossover from 
the power-law correlations of the 5=1/2 chain to the 
long-range order of the square lattice. Interestingly, lad- 
ders with an odd number n c of chains have power-law 
spin correlations in their ground state, while those with 
even n c exhibit exponentially decaying correlations due 
to the presence of a spin gap || . In analogy to the general 
5 chain [Q] , the fundamental difference between even and 
odd ladders is thought to be of topological nature |||| ■ 

In this Letter we investigate isotropically coupled lad- 
ders as well as the two-chain ladder in the regime of 
weak inter-chain couplings. The former systems are di- 
rectly relevant to the known experimental systems which 
are approximately isotropic p|. Moreover, the crossover 
from ID to 2D is most naturally studied in this case. 
We perform numerical simulations which allow us to de- 
termine both the correlation length and the static struc- 
ture factor down to very low temperatures. Our results 
provide a basis for comparison with future neutron scat- 
tering and NMR experiments. The two-chain ladder at 



weak inter-chain couplings allows us to investigate the 
formation of a spin gap away from the unstable T = 
fixed point of the gapless 5 = 1/2 chain [[To), and in- 
deed serves to demonstrate the extreme fragility of the 
5 = 1/2 chain power-law correlations. Our data for the 
intra- and inter-chain correlation functions at low tem- 
perature arc in excellent agreement with a recent theo- 
retical prediction JO]. Moreover, we discover that in this 
regime the correlation length exhibits universal scaling 
behavior with respect to the inter-chain coupling. 

The Hamilton operator for a Heisenberg system of n c 
5 = 1/2 chains forming a ladder is 



(1) 



Here, S, 



is the quantum spin operator located at 



the point i, while and (ij)± denote nearest neigh- 

bors along and between chains, respectively. We consider 
antiferromagnetic couplings, that is Jn, J± > 0, and peri- 
odic boundary conditions along the chains. We use units 
in which Ti = ks — 1 and, unless noted otherwise, J|| = 1. 

The ladder systems are investigated with a very effi- 
cient loop cluster algorithm [|l2|Jl3| 1 which allows access 
to very low temperatures and the implementation of im- 
proved estimators in order to reduce statistical errors 
of observables. We simulate lattices large enough, both 
along the chains and in Euclidean time, so that finite size 
and finite Trotter number effects are comparable to the 
statistical errors. In particular, the length of the ladders 
is kept ~ 15 times larger than the correlation length. 
Typically, 4 x 10 4 loop updates are performed for equili- 
bration, followed by 4 x 10 5 measurements. 

In order to obtain information about the gap A(n c , J±) 
for even n c , we measure the uniform susceptibility 

X (n c ,J ± ;T)~T- 1 ((^5f) 2 ). (2) 

i 

Gap values extracted from fits to the low-T form Jill 



1 



x(n c ,J±;T)~r- 1 /2 e -A(n.,Jx)/T 



(3) 



are shown in Fig. 1. For J± > 2, we find very good 
agreement with strong-coupling expansion results for 
n c = 2 and 4 Ea], shown as solid lines. Our suscepti- 
bility data for the isotropically coupled ladders (Jj_ = 1) 
agree well with recent Monte Carlo work [M, and we 
obtain A(n c ,l) = 0.502(5), 0.160(5), and 0.055(6) for 
n c = 2, 4, and 6, respectively. In an earlier study 
@, A(2,l) = 0.504 and A(4, 1) = 0.190 were found. 
For n c = 2, we are able to access the weak-coupling 
regime characterized by A (2, Jj_) ~ J± 0,0. We find 
A(2, J±)/J± = 0.41(1), which is somewhat smaller than 
A(2, J±)/J± = 0.47(1) obtained previously p7 . 
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FIG. 1. Dependence of the spin gap on the inter-chain 
coupling. The solid lines represent strong-coupling results 
[Jl5| , and the dashed line indicates the weak-coupling behavior 
A(2, J x -f 0) = 0.41(1) Jx. 

Next we compute the staggered correlation function 



C(i,j) = (-l) sl S n ^)(S,-S i ), 



(4) 



where sign(i, j) is 1(— 1) if the spins at i and j are sepa- 
rated by an even (odd) number of couplings. For a system 
of two weakly coupled chains at T = 0, conformal field 
theory predicts jllj 

C(r) = G+ (r)G_ (r) [G_ (r)GL (3r) ± G+ (r)G+ (3r)] , (5) 

and this form is claimed to be exact in the continuum 
limit at small J±. The sign of the last term is plus for 
intra- and minus for inter-chain correlations, and r mea- 
sures the distance of two spins along the ladder. The 
functions 



G±(r) = r~ 1/4 F ± (r/0[l ± 2~ 3 / 2 C 1 } + 0(r 



-5/4 



(6) 



are correlation functions of the 2D Ising model, with scal- 
ing functions F± |Q. In Eq. (6) £ denotes the spin-spin 
correlation length of the ladder system. The intra- and 



inter-chain correlation functions differ only at short dis- 
tances. At large distances r/£ 3> 1, C(r) decays as 



C(r) 



-r/£ 



(7) 



with A = 1/2, which is equivalent to the 2D Ornstein- 
Zernike (OZ) form. Figure 2 shows the low-T correlation 
function for n c — 2 at a weak coupling of Jj_ = 0.2. The 
lines are the result of a fit to Eq. (5) with only two 
fitting parameters: £ and an overall amplitude. The fit is 
excellent over the entire range, and we obtain £ = 19.4(2). 
As is evident in Fig. 2, Eq. (5) correctly captures the 
crossover from the short- to the long-distance behavior 
with a concomitant change in length scales by ~ 3. 
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FIG. 2. Intra- and inter-chain correlation functions at 
T = 0.01 for a, S — 1/2 two-chain Heisenberg ladder of length 
300 with Jx = 0.2. Note that periodic boundary conditions 
were employed along the chains and that for clarity C(r) is 
shown only for even r. The lines are the result of a fit to Eq. 
(5) in the symmetrized form [C(r) + C(300 — r)]. 

In order to deduce £(n c , Jx; T) for general n c , Jj_, and 
T, we fit C(i,j) at large distances r > 3£ to the asymp- 
totic form Eq. (7). For n c = 4 and 6, Eq. (7) with 
A = 0.5 describes our low-T data very well. However, 
we find that the asymptotic behavior for even n c crosses 
over between T w 0.2A(n c , Jj_) and T « 0.4A(n c , J±) to 
the ID OZ form, that is, A = 0. For odd n c we find that 
the ID form works very well at all temperatures. 

In Fig. 3, the result of this analysis for J± = 1 is 
shown together with £(T) for the square lattice as ob- 
tained by both Monte Carlo |L9| and neutron scatter- 
ing in Sr-iCuOiCl2 Due to the presence of a gap 
for even n c , £(T) remains finite in the limit T — > 0. 
We estimate that £(n c ,l;0) = 3.24(5) and 10.3(1) for 
n c = 2 and 4, respectively. The result for n c = 2 is 
very close to the value £ = 3.19(1) obtained previously 
H. However, in Ref. || £ = 5 to 6 was obtained for 
n c = 4, which is significantly smaller than our result. 
For n c = 2 and 4 we obtain the respective velocities 
c(n c , 1) = A(n c , l)£(n c , 1; 0) = 1.63(2) and 1.65(3) which 
lie in between the ID and 2D values c\d = 7r/2 and 



2 



c 2D ^ 1-68. 

The correlation length of the single S — 1/2 chain has 
been determined in a thermal Bethe-ansatz study ]2Q|: 



Cid(T) ~ T [2 - (1 - 0.4866" 1 ln(6))] 



(8) 



with b = -ln(0.3733T). For T < 0.3 our data for the 
chain (n c = 1) agree with this low-T form, indicated 
by the dashed line in Fig. 3. We observe that with 
decreasing temperature £ for n c = 3 and 5 gradually 
approaches 
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FIG. 3. Correlation length of n c isotropically coupled anti- 
ferromagnetic S = 1/2 chains. The result for the single chain 
(n c = 1) is compared with the theoretical prediction Eq. (8) 
(dashed line). Also shown are the results for the square lat- 
tice obtained from both Monte Carlo simulations Jl9| as well 
as neutron scattering experiments in Sr^CuO^Cl^ |H. 

In Fig. 4 the static structure factor at (7r,7r), C n>ir = 
is shown for the isotropically coupled sys- 
tems. The overall trend with n c is similar to that for £. 
However, C 7ri7r (n c , Jj_;T) for even n c exhibits a peak at 
a temperature T max well below the corresponding spin 
gap. Not surprisingly, T max appears to coincide with the 
temperature at which the asymptotic correlation func- 
tion Eq. (7) begins to cross over from the ID OZ form 
at higher temperatures to the low-T 2D form. 

Figure 5 shows the correlation length for the two-chain 
ladder at several inter-chain couplings J± < 1. For 
T w A(2, Jj_), f(2, Jx;T) is larger than£i£)(T), since the 
ladder system has a higher effective coordination number 
than the chain. Well below T = A(2, Jj_), the pres- 
ence of the non-zero spin gap becomes apparent, and 
£(2, J±;T — > 0) remains finite. With decreasing J±, the 
gap decreases as shown in Fig. 1, and correspondingly, 
£(2, J x ;0) increases from 3.24(5) at J x = 1 to 25.6(2) 



at J± — 0.15. In the weak-coupling regime, Jx < 0.5, 
the velocity c(2, Jx) = A(2, Jx)£(2, J±; 0) equals that of 
the single chain, cirj = 7r/2, to within the error of our 
analysis. In this regime the heuristic form 



£-'(2, Jx; T) = cri A(2, Jx) + &(T)e 



-A(2,J±)/T 



(9) 



describes our correlation length data very well, as shown 
by the solid lines in Fig. 5. Note that since the spin gap 
has been determined from measurements of x(T), this 
comparison contains no free parameters. 
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FIG. 4. Static structure factor at (n, n) for isotropically 
coupled antiferromagnetic S =1/2 ladders, as well as for the 
single chain (n c = 1) and the square lattice |Ti|. For even n c , 
Cr.ir exhibits a peak (indicated by arrows) at a temperature 
below that corresponding to the relevant spin gap. 

The temperature dependence of £(2,Jx;T) in the 
weak-coupling regime primarily results from that of the 
single chain, Eq. (8). Apart from logarithmic correc- 
tions, the latter is simply £ld(T) ~ T~h Eq. (9) thus 
becomes (A^)- 1 = 2/7r + ^4(T/A)e _A/T , which suggests 
plotting A£ versus T/A to test for the anticipated scal- 
ing for Jx < 0.5. As shown in the inset of Fig. (5), our 
correlation length data indeed collapse onto a universal 
curve. Over the indicated range the effective value for A 
is ~ 1.7 compared to the low-T value of 2. 

For the two-chain ladder, we have established that 
A(2, Jx — » 0) = 0.41(1) Jx- We can furthermore estimate 
the weak-coupling behavior for n c = 4. From Fig. 1, it 
can be inferred that A(4, J±) ~ A(2, Jx) — 0.35Jx, which 
leads to A(4, J± — > 0) = 0.06 Jx as a lower bound. From 
the lowest J ± for n c = 4, we have A(4, 0.6) = 0.12(1) Jx 
which constitutes an upper bound. Thus we arrive at the 
estimate A(4, J ± -> 0) = 0.09(3) Jx- 

We note that the trend with increasing n c of the ladder 
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gap resembles that of the single integer-spin chain as a 
function of increasing 5: For S — 1, Ag = i = 0.4105J|| 
is known very accurately |21|, while recent estimates for 
5 = 2 range between A s=2 = 0.049(18) J\\ and 0.085(5) J\\ 
PH . It has recently been argued that the two-chain lad- 
der with antiferromagnetic inter-chain coupling lies in 
the same phase as the 5=1 single chain p3| , and that 
the weak-coupling physics is independent of the sign of 
J± jll]]. Clearly, further numerical and analytical work 
to establish this correspondence is necessary. 



of the correlation length are in excellent agreement with 
Monte Carlo simulations and theory , as well as high- 
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FIG. 5. Correlation length of two antiferromagnetically 
coupled S = 1/2 chains for J± < 1. The dashed line is 
Eq. (8), and the solid lines are Eq. (9) with A(2, J±) as 
obtained from measurements of x{T)- Inset: Scaling plot of 
A(2, Jl)£(2, Jj_;T) versus T/A(2, J ± ) for J x < 0.5. 

In summary, our Monte Carlo calculations on spin lad- 
ders have provided a number of new results. First, our 
measured low-tcmpcrature correlation function for the 
weakly-coupled two-chain system agrees very well with 
the predicted theoretical form fllj, Eq. (5), which is 
claimed to be exact. In particular, Eq. (5) captures 
the crossover in behavior between short and long dis- 
tances. With increasing temperature the even ladder cor- 
relations cross over from the 2D Ornstein-Zernike form 
to the ID form of pure exponential decay. Second, for 
J± < 0.5 we have discovered that A(2, Jj_)£(2, J±;T) 
exhibits a rather simple scaling behavior in the variable 
T/A(2, Jj_). Third, the structure factor for even n c ex- 
hibits a peak at a temperature well below the correspond- 
ing spin gap. Our results for the correlation length and 
the static structure factor of the isotropically coupled 
ladders pertain to the experimental systems. For the 
prototype 5 = 1/2 square-lattice Heisenberg antiferro- 
magnet Sr2CuC>2Cl2, neutron scattering measurements 



temperature series expansion 
crystals of the ladder systems become available, we in- 
tend to extend our neutron scattering work to these in- 
teresting systems. We also hope that our Monte Carlo 
results will motivate future theoretical efforts to predict 
C(i,j) for general n c , J±, and T. 
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